Effective screening of localized charged perturbations 
in metallic nanotubes: roles of massive bands 
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The massive-band effects on screening beliavior of metallic carbon nanotubes are theoretically 
investigated using two different methods; continuous and lattice quantum theories. Both approaches 
show screening of a localized external perturbation with an effective screening length of the order 
of the nanotube diameter. Calculating the nonlinear deformation of the local density of states near 
the charged perturbation, we show that the perturbative effects of the massive bands are effectively 
canceled by direct massive band interactions, such that a good agreement between the two methods 
can be achieved. The effective screening is important in nanoscale integration of nanotube-based 
electronic devices. 
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I. INTRODUCTION 

The response of the electrons in carbon nanotubes 
(CNT's^'^) to external charged perturbations is governed 
by the Coulomb interactio n'^i"^'^ . Several aspects of the 
Coulomb interaction are already seen in some experi- 
ments of chargingSiS and of temperature dependent re- 



sistivity- which are consistent with current understand- 
ing of low-energy excitations in metallic CNT's. The 
low-energy excitations in the metallic CNT's are theo- 
retically modeled by massless bands (linear dispersion 
bands) with the Coulomb interactionSi^. However when 
we consider very localized perturbation, it is not clear if 
the low energy theories describe the physics correctly. 
Because short distance physics corresponds to high en- 
ergy, massive bands (other sub-bands except the massless 
bands) would come into dynamics From the applica- 
tion point of view, understanding the response of CNT's 
to localized perturbations (e.g., individual dopant atoms) 
is crucial as it provides the extents of the depletion region 
in nanotube-based electronic devices. The importance of 
this issue is observed by considering the fact that in any 
nano-electronics application of nanotubes, different func- 
tional parts of the device, e.g., drain, source, etc., should 
be doped differently. In order to achieve nanoscale in- 
tegration, these differently doped parts should be sepa- 
rated from each other by distances of the order of a few 
nanometers. This would not be possible if the depletion- 
layer length is not of the order of a few nanometers itself. 
In this letter, taking massive band effects explicitly into 
account, we focus on the screening of an external charge 
localized on the surface of the metallic CNT's. 

It is well-known that in lower dimensional materials 
in which the internal electrons can move only in one 
or two dimensions, asymptotic behavior of the poten- 
tial produced by an external charge does not become an 
exponential as a function of distance from the external 
charge. Examples of these lower dimensional materials 
are graphite sheet and carbon nanotubes. DiVincenzo 
and Meleii analyzed the induced charge around an ex- 
ternal charge on the graphite sheet and found that the 



external charge attracts some induced charge with a char- 
acteristic length scale of the system which they called "ef- 
fective screening length" . A similar definition is used to 
describe the effective screening in metallic carbon nan- 
otubes by Lin and Chuu^". In the present paper, we 
study the details of the pattern of induced charge around 
an external charge on the metallic nanotubes in order to 
examine the effects of massive bands on screening. 



The high-energy excitations can in principle affect the 
screening phenomena through two different mechanisms; 
adding massive-fermion loop corrections to the interac- 
tions of massless bands and direct interaction of massive 
bands. Two methods to deal with these mechanisms are 
investigated. One is a quantum field theoretical method 
dealing with the massless bands and the one-dimensional 
long-range Coulomb interaction^^, with one- loop vacuum 
polarization correction by the massive bands. The other 
is a lattice model employing a fully self-consistent tight- 
binding approach including the interactions of massive 
bands^'^. It is the purpose of the present work to compare 
the results of these two different methods and investigate 
the screening nature in CNT's. Especially we shall com- 
pare the induced charge distribution around the external 
charge. Our results are advantageous compared to the 
previous works on screening in metallic nanotubesiS, as 
we include an effective thickness for the nanotube wall, 
and derive the dependence of our results on the corre- 
sponding cut-off parameter. This inclusion is crucial in 
reaching an agreement between the continuous and lat- 
tice results. Additionally, in our self-consistent tight- 
binding model there is no linear response assumption, 
and the deformation of local density of states (LDOS), 
due to the creation of localized states, are included. We 
also explicitly show how high-energy electron/hole exci- 
tations respond differently to the presence of a charged 
perturbation. 
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II. CONTINUOUS FIELD THEORY: ONE-LOOP 
MASSIVE FERMION CORRECTIONS 



We begin with the quantum (continuous) field theo- 
retical approach to the screening eff'ect. The 7r-elcctrons 
in CNT's are described by two-component fermion fields 
ipNs specified by the band index N{£ {1, • ■ •}) and spin 
{TjI})- The Hamiltonian density of this system is 
given hy H — Hf + Ti-c where is the free (unper- 
turbed) part and Tic is the perturbed part including the 
Coulomb interactions. If we include only tight-binding 
interactions, then Hp = '^ns'^n where 
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is the hopping integral, ky = —idy is the wave vector 
along the tubule axis, the wave vector around the circum- 
ference ak^ = — -2^ -I- is discretized by the periodic 
boundary condition where d is the diameter of the CNT 
and a the nearest-neighbor C-C spacing. Here we take 
the case of metallic zigzag CNT's in which the maximum 
number of N is A^o (number of hexagons along the cir- 
cumference); armchair case is analyzed similarly. The 
function z{kx, ky) is given by z{kx,ky) = ^ ■ e''^'"' . The 
vector Ui is a triad of vectors pointing in the direction of 
the nearest neighbors of a carbon site^^ . The free Hamil- 
tonian densities Hn^/^ and 'H^Nq/s describe massless 
fermions (electrons in massless bands) and the others are 
for massive fermions in the metallic zigzag CNT. The en- 
ergy eigenvalues of Hat are given by ±t4|z(A:^, /cy)| which 
reproduce a standard band structure^. These fermions 
interact with each other via the long-range Coulomb in- 
teraction: 



J{y)V{y - y')J{y')dydy\ 



(2) 
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where V{y) is the unscreened Coulomb potential that 
we will specify shortly, and the internal charge density 
J is given by the sum of individual fermion densities: 
J = X^ATs "^Ns where Jjvs = ipl^^ip^s- Here we denote the 
integral region D = [0 : L] and L is the system length. 

Low-energy excitations in metallic CNT's are de- 
scribed by a field theory of four massless fermions in- 
teracting with each other and with massive fermions 
through the long-range Coulomb interaction. The mas- 
sive fermions modify the Coulomb potential between 
massless fermions because of the vacuum polarization. 
We estimate this effect using one-loop perturbation. 
The modified Fourier component of the Coulomb po- 
tential is given by /?„ — , where /?„ is the n-th 

Fourier component of the unscreened Coulomb potential, 
V{x) = E„/^ne"'2™x/L^ j^j^^ Tr"(< 0) in the denom- 
inator denotes the one-loop contribution of the massive 
fermionsi^. The unscreened Coulomb potential along the 



Vix) = 
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with K{z) being the complete elliptic integral of the 
first kind and Oz = 1.3 [A] a cutoff length which cor- 
responds to the Hubbard coefficient Uh = ^V-r: in the 
lattice model (Eq|SJ). The massive fermions are ob- 
served to weaken the Coulomb interaction among mass- 
less fermions as is naturally expected. Because the 
Coulomb potential is modified by the vacuum polariza- 
tion effects, the massless fermions interact with each 
other via the screened Coulomb potential whose Fourier 
components are given by /?„. 

An external charge can be included by replacing J with 
J -f J™ in the Coulomb interactions (Eq. [SJl, where J'^^ 
is a c-number whose integral over the nanotube length 
is equal to the total external charge. Here we take a 
point external charge distribution which is modeled by 
J'^^ = 5{x — xq) where xq is the position of the external 
charge. To calculate the induced charge density, we use 
one-loop approximation which gives an exact result in the 
case of strictly linear dispersion relationi^. In Fig. the 
induced charge distributions around the external point 
charge for both the unscreened potential and the screened 
one are depicted. These are obtained using the following 



n>0 



2 /27m, ^, 
cos ( ^— (x - xo) , (4) 
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where T™' is the one-loop amplitude of massless bands. 
It should be noted that the effective screening length is 
about the diameter of the nanotube for both cases and 
is not sensitive to the system lengtbi^. 

What we ignore in the above formula of induced charge 
density is (a) higher-order perturbative effects of mas- 
sive loops on the Coulomb interaction between massless 
fermions, and (b) induced charge due to massive fermions 
(direct massive-band effect). Lin and Chui*iS included 
the latter effect by adding T^' to T^' in the numera- 
tor of Eq.Q. This assumption, however, can only re- 
sult in a linear relationship between the external and 
induced charges, while it is shown that linear-response 
theory is not sufficient at least in the case of screening in 
graphiteii. 



III. SELF-CONSISTENT TIGHT-BINDING: 
NONLINEAR SCREENING 

Up to now, only the massive-fermion corrections to the 
massless-hand interactions, through the leading terms of 
vacuum polarization, were considered. We now turn to 
the massive-hand effects on the screening phenomena, 
and investigate the inffuence of direct massive-band inter- 
actions. This is accomplished by using a lattice model to 
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FIG. 1; Total induced charge (including spin degeneracy) 
around the point external charge for the (5,5) armchair tube. 
The dotted line is obtained from the lattice model. The 
solid and dashed lines are the continuous field theoretical re- 
sults. The dashed one is for the calculation using unscreened 
Coulomb potential /3„ , and the solid one is for • The screen- 
ing length is about the diameter of the nanotube. Similar 
results are obtained for the zigzag (9,0) tube. 

describe the nanotube. The model is essentially a tight- 
binding onei^ which is generalized by including long- 
range Coulomb interactionsi^. 

Within this model, the Hamiltonian of a perfect, infi- 
nite nanotube is written as H = Hp + He with 

Hf = Y,V.c\c,, (5) 

Hc = Y. //"^ _, 4c. (6) 

where Uh is the Hubbard parameter used in describing 
on-site interactions. Using the experimental estimation 
oiUn for carbonic, we set Uh = 4T4- Uh (corresponding 
to to the cut-off flz) is the only parameter of our model 
and we examine the dependence of our results on this 
parameter, by considering several cases of Uh (or a^). In 
EqEl Sui = Hi — is the change in the self-consistent 
and screened occupation number Ui at site i, as compared 
to the occupation imposed by external perturbation, n°. 
Therefore, Srii would be the induced charge at site z, in 
response to the external perturbation. We are concerned 
with localized perturbation, which is modeled by modi- 
fying n° as n° — > n° + . 

We distinguish three different parts in the infinite nan- 
otube: An unperturbed semi-infinite part to the left, a 
perturbed finite part in the middle, and an unperturbed 
semi-infinite part to the right. In the unperturbed parts, 
— 0.5 (excluding the spin degeneracy) and Srii — 0. 
This implies that the length of the perturbed finite part, 
for which Jrii's are assumed to be nonzero, is chosen to be 
long enough such that the external perturbation is fully 
screened within this part of the nanotube. For concrete- 
ness, we assume that the localized external perturbation 



FIG. 2: The cutoff, Uh (or a^), dependence of the ratio de- 
fined by D/W , where D is the peak of the induced charge 
at the position of the external charge and W is the one-half 
of the width at the half peak (see Fig. . The results show 
that bigger Uh (smaller a^) gives smaller screening length. 
Experimentally, Uh was reported-- to be between 2.4 and 4 

only modifies the sites of the two carbon rings of the nan- 
otube by adding nf^'s to the corresponding n°'s. (A ring 
consists of all the carbon atoms with the same longitudi- 
nal coordinate along the nanotube axis.) The two rings 
with nonzero nf^ 's are assumed to be at the middle of the 
perturbed finite part of the nanotube, which, in turn, is 
taken to consist of 16 carbon rings. This number of car- 
bon rings has proven to be long enough in order to 
derive the converged induced charges, as increasing the 
number of rings has had no effect on the self-consistent 
occupations. 

The self-consistent procedure of calculating the in- 
duced charges (Jn^'s proceeds as follows ~: Starting with 
an initial guess of zero Srii's in the Hamiltonian Eq|Sl 
the surface Green's functions, Gl and Gr, of the two 
semi-infinite unperturbed parts to the left and right of 
the perturbed region, are calculatedi^. These surface 
Green's functions are then attached to the Green's func- 
tion of the perturbed middle part, Gp = {z — Hp)~^, 
where z is the complex energy and Hp is the Hamilto- 
nian of the perturbed region, in order to obtain the total 
Green's function of the system, Gt, projected onto the 
perturbed regioni2i22i2iiS2i2i^: 

Gt - (Gpi - VplGlVlp - VprGrVrp)-\ (7) 

Here, V matrices indicate the coupling between the per- 
turbed region and the first unit cells of the unperturbed 
parts to the left and right. Using Eq[7| the LDOS at site i 
in the perturbed region, gi{E) = [Grih i) — G'fp{i, i)]/2iT, 
is derived. The induced charges can be obtained through 
integrating gi[Eys. Srii^s are then used to initialize the 
Hamiltonian for another round of Srii's calculation. This 
procedure is continued until the induced charges are de- 
termined self-consistency. Throughout the self-consistent 
calculation, the chemical potential of the system is as- 
sumed to be that of the unperturbed nanotube, as after 
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screening the whole system would be in thermal equilib- 
rium and would have a unique chemical potential. Fur- 
thermore, the external perturbation being localized and 
finite, the unique chemical potential would be that of the 
unperturbed semi-infinite parts of the nanotube, which 
act like reservoirs attached to the perturbed finite part. 

As specific examples, an armchair (5,5) and a zigzag 
(9,0) nanotube are considered, whose diameters are 6.78 
and 7.04 [A]. The external perturbation of the middle two 
rings of the perturbed part of these tubes are defined as 
n^^ = 0.0250 and nf^ — 0.0278, for the armchair and 
zigzag tubes, respectively. This external perturbation 
amounts to adding a total charge of l.Oe to the middle 
two rings for both the armchair and the zigzag tubes. 

Comparing the self-consistent procedure used in the 
tight-binding model of the present work with the tight- 
binding approach of Lin and Chuijiii, we notice that in 
our model there is no restriction on the shape of the 
LDOS within the perturbed region. In fact, the very 
shape of the LDOS is calculated self-consistently, as de- 
scribed above. Lin and Chun's approach, however, as- 
sumes fixed energy dispersion relations for the whole nan- 
otube, and restricts the modifications due to the charged 
perturbation to the linear response theory via dielectric 
constant and response function. Our model is therefore 
advantageous, as it properly takes into account the non- 
linear effects of the localized states which are produced 
within the perturbed region. This is explicitly observed 
from the LDOS curves in the perturbed region, as will 
be shown shortly. 



IV. COMPARISON AND DISCUSSION 

The resultant induced charges for the (5,5) armchair 
tube are depicted in Fig. It is observed that the ef- 
fective screening length is of the order of the nanotube 
diameter, in all the different calculations. The length also 
depends on the cutoff or the Hubbard coefhcient Uh- 
The cutoff dependence of induced charge pattern is plot- 
ted for several Uh values in Fig. [51 The main features 
of the induced charge pattern are observed to hold for 
the whole range of experimentally reported Uh^- From 
Fig. n] we further see that the induced charges obtained 
by including only the massless bands in the continuous 
field-theoretical calculation are in qualitative agreement 
with the results of lattice calculation, which, in addi- 
tion to the massless bands, includes the effects of direct 
massive bands interactions. Interestingly, including the 
massive-loop corrections to massless bands interactions 
does not result in better agreement with the results of 
lattice model. This suggests that, within the continu- 
ous model, the massive-loop corrections to the massless 
band interactions somehow cancel the effects of direct 
massive band interactions, such that the results of a cal- 
culation based only on massless bands turn out to be 
in agreement with the results of the lattice model. In 
other words, in the continuous approach, the induced 
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FIG. 3: Self-consistent electron (top) and hole (bottom) 
LDOS. Fermi energy is located at zero. 



charge density is given by the sum of two densities as 

(J) = (Jmasslcss) + ( Jmassivc) whcrC (Jmassless) is defined 

in Eq^ (^massive) is an induced charge density due to 
massive bands themselves, and is what we have ignored. 
Perturbative inclusion of massive bands reduces the peak 
of the former charge density and results in a more ap- 
preciable difference with the lattice results. Considering 
overall agreement of the lattice results and continuous 
field theoretical results based on the massless bands, it 
is natural to expect that (Jmassivo) partially cancels the 
perturbative effects of the massive bands on the mass- 
less bands interaction, such that the net induced charge 
density is basically decided by the massless bands. 

In addition to the calculations for the armchair tube, 
we have performed similar calculations on the (9,0) zigzag 
tube. The results of the continuous modeli2. show the 
same features, as those of the armchair one, when com- 
pared to the results of lattice model. This indicates that 
the agreement between the continuous massless results 
and lattice results is not accidental. It is remarkable that 
in both the continuous and lattice calculations the peak 
value of the induced charge shows a reduction of ~ 15%, 
when one compares zigzag results with armchair results. 

In order to see the LDOS deformation, the creation 
of localized states within the perturbed region, and the 
effects of massive bands more clearly, we show the self- 
consistent LDOS's for a typical unperturbed site (solid 
line), a site located at one of the rings perturbed exter- 
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nally (dashed line), and a site, of the perturbed portion 
of the nanotube, located at the furthest ring relative to 
the middle two rings (dotted Une) in Fig. El The lat- 
ter site indicates the border between the perturbed finite 
part of the tube and one of the unperturbed semi-infinite 
parts. In Fig. enlarged typical regions of the electron 
(hole) LDOS's are depicted. It is clearly seen that the 
densities of massive states in the perturbed region are 
different from the corresponding densities of the unper- 
turbed region. However, as one expects, the difference 
is much more pronounced near the location of external 
perturbation than far from it. This is indeed an indica- 
tion of localized-states creation. The border LDOS shows 
how the extremely perturbed LDOS near the external 
perturbation transforms rather smoothly, as one moves 
toward the border of the perturbed region, into the un- 
perturbed LDOS. Another very interesting feature is the 
difference in the behavior of the electron states as com- 
pared to the hole states. For the external perturbation 
mentioned above, namely addition of one extra electron 
to the middle two rings, it is seen that the hole states 
tend to accumulate at lower energies, while the electron 
states are pushed toward higher energies. This indicates 
the decisive effect of the Coulomb interactions among the 
massive fermion- (both electron- and hole- ) states. One 
should notice that, due to the rotational symmetry selec- 
tion rule the electron/hole states belonging to a cer- 
tain band are not allowed to jump over to another band 
with higher/lower energy. But within each band the elec- 
tron states tend to concentrate at higher energies, and the 
hole states tend to concentrate at lower energies. This is 
indeed the reason for the enhanced/decreased densities 
at band edges. 



V. CONCLUSIONS 

In summary, we have analyzed the charge screening 
in metallic nanotubes using two different methods, 
including explicitly the massive bands effects, based 
on lattice and continuous models. In both methods, a 
localized external charge is shown to be screened with 
an effective screening length about the diameter of the 
tube. The local density of states, calculated by the 
lattice model, indicates the formation of the localized 
states within the perturbed region, and shows that the 
massive bands are strongly affected by the external 
charge. As for the continuous field theoretical model, we 
have included the massive bands one-loop corrections in 
massless bands interactions. Although non-perturbative 
effects of massive bands were not considered, the 
results show basically the same pattern as that of the 
lattice model results. This shows that the perturbative 
corrections of massive bands are efficiently canceled by 
the induced charge due to massive bands themselves. 
The short effective screening length reported here makes 
it possible, at least theoretically, to achieve nanoscale 
integration of nanotube-based electronic devices. 
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